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Abstract 


We consider the case of a scalar field, the inflaton, coupled to both lighter 
scalars and fermions, and the study the relaxation of the inflaton via particle 
production in both the linear and non-linear regimes. This has an imme¬ 
diate application to the reheating problem in inflationary universe models. 
The linear regime analysis offers a rationale for the standard approach to the 
reheating problem, but we make a distinction between relaxation and ther- 
malization. We find that particle production when the inflaton starts in the 
non-linear region is typically a far more efficient way of transfering energy out 
of the inflaton zero mode and into the quanta of the lighter scalar than single 
particle decay. For the non-linear regime we take into account self-consistently 
the evolution of the expectation value of the inflaton field coupled to the evo¬ 
lution of the quantum fluctuations. An exhaustive numerical analysis reveals 
that the distribution of produced particles is far from thermal and the effect 
of open channels. In the fermionic case, Pauli blocking begins to hinder the 
transfer of energy into the fermion modes very early on in the evolution of 
the inflaton. We examine the implications of our results to the question of 
how to calculate the reheating temperature of the universe after inflation. 
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I. INTRODUCTION 


Thermalization, reheating and relaxation towards equilibrium are ubiquitous non¬ 
equilibrium phenomena that play a very important role in the early universe, in heavy 
ion collisions and in the dynamics of phase transitions away from equilibrium. 

In particular in typical inflationary scenarios of the early universe, after the many e- 
folds of inflation necessary to solve the horizon and homogeneity problems, the matter 
and radiation energy density had been red-shifted to almost zero. However, once inflation 
has ended, this scenario has to merge with the standard big bang, radiation dominated 
cosmology. The succesful intertwining of inflation with the standard big bang cosmology 
thus necessitates some source of energy and entropy to rethermalize, or reheat the universe 


The standard picture of old and new inflation invokes a scalar held, the inhaton, that 
produces a phase transition (either second or first order) and that at the end of inflation the 
expectation value of its zero momentum mode oscillates about the minimum of its (effective) 
potential. In chaotic inflation ||, the inhaton energy density still drives an inhationary 
phase, but there need not be any phase transition for this to occur. This inhaton held is 
constrained |I|-[3| to have a mass a few orders of magnitude smaller than the Planck mass, 
but typically much larger than the masses of the particles involved in the particle physics 
models. In the standard view of inhationary models, the expectation value of the inhaton 
held oscillates around the minimum of its potential, and its couplings to the lighter particles 
allow the inhaton held to decay. This decay process is then supposed to induce a damping 
term in the evolution equation for the inhaton expectation value of the form T</> with T being 
the total decay rate of the inhaton held BS# The standard estimate for the reheating 
temperature based on single-particle decay @ is then obtained by comparing the total decay 
rate of the particle (T) to that of the expansion, obtaining T r ~ 0.1 ^TM p i |2|]. 

Recent investigations of the non-linear quantum dynamics of scalar helds reveal a variety 
of new and striking phenomena [Q TA\. The main relevant implication for the reheating 
problem being that the particle production induced by the time evolution of the inhaton is 
signihcatively different from linear estimates. 

The non-linear (quantum) effects lead to an extremely effective dissipational dynam¬ 
ics and particle production even in the simplest self-interacting scalar field theory [Til.[hill hi 
which single particle decay is kinematically forbidden. In the case in which only the classical 
evolution of the expectation value is taken into account in the evolution of the quantum fluc¬ 
tuations, this mechanism corresponds to parametric amplification or parametric resonance, 
because the classical time evolution is periodic in time |T5|. However when back-reaction 
effects are taken into account, particle production damps out the evolution of the scalar 
field, and this damped evolution is incorporated in the equations for the quantum fluctua¬ 
tions. In this case the time evolution of the scalar field including the dissipative effects is not 
periodic; this implies that the situation is far different from parametric amplification. We 
refer to this situation, which includes the back-reaction effects as induced amplification, 
to distinguish it from parametric resonant amplification [ |T5| . 

Back-reaction effects |ll],T3|] drastically change the picture of particle production. In the 
case of parametric amplification, particle production never shuts off and the total number of 
particles created is infinite in Minkowski. It is only when the back-reaction is incorporated 
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that the damping effects on the expectation value feed back to the particle production 
mechanism, eventually shutting it off. This has been studied in detail analytically and 
numerically in [|T3| for a self-interacting scalar field. 

Induced amplification is most effective when the amplitude of the expectation value of the 
scalar field (from the equilibrium value) is large and is, therefore, an intrinsically non-linear 
process, which we refer to as non-linear relaxation. 

Particle production via parametric amplification of quantum fluctuations has been stud¬ 
ied recently in the semiclassical approximation, but without back-reaction effects, in con¬ 
nection with fermion production by (pseudo) Nambu-Goldstone bosons [TTh . 


Induced amplification and non-linear relaxation will be the primary mechanism for dis¬ 
sipation via particle production in the case of large amplitude of the scalar field. Such is 
the case in chaotic infation scenarios [5|, as well as in the out of equilibrium regime during 
phase transitions for fields with very flat potentials near the maximum, which is typically 
required for a long inflationary phase and also occurs for the moduli fields in string theory 

Q 

As anticipated by Kofman, Linde and Starobinsky |Ll|, this mechanism will necessarily 
modify the standard picture of reheating that was primarily based on the premise that 
particle production only ocurred when the inflaton oscillates with small amplitude at the 
bottom of the (effective) potential. Understanding the time scales and non-linear processes 
of reheating and eventually thermalization acquires further importance with the possibility 
that asymptotic oscillations of the inflaton field around the minima of the (effective) potential 
may still be present in today’s universe in the form of dark matter |TT1,[T8|] . 

Thermalization is a process that is fundamentally different from particle production. 
Typically, as will be seen below, the particles produced by the process of induced amplifi¬ 
cation will be in a non-thermal distribution [[□]]. Thermalization is a collisional process, in 
which the produced particles will exchange energy (and momentum) and eventually achieve 
a thermal distribution. In principle the time scales for dissipation via particle production 
may be widely different from the time scales of collisional thermalization. For very weakly 
coupled theories, if the distribution of produced particles via induced amplification is very 
far from thermal, many collisions will be necessary to thermalize the system and the time 
scale for thermalization may very well be large compared to the time scale of dissipation via 
particle production. This will have implications in cosmology discussed in the conclusions. 

The process of thermalization, reheating and relaxation of perturbations, and the ensuing 
production of entropy is also very relevant in heavy ion collisions |T^| and in phase transitions 
in particle physics, at the electroweak scale within the context of baryogenesis [2C)j 
as for the quark-gluon plasma, deconfining and chiral phase transitions 


as well 


Our main goal in this article is to study the time scales associated with the process of 
non-linear relaxation and compare these to those of single particle decay which we refer to 
as “linear relaxation”. 

In the linear relaxation case the important scale is determined by the decay rate T, usually 
referred to as the “damping rate”. The literature usually identifies T sa ImE(w,fc)/w 
with the “thermalization rate”. Here ImE(w, k ) is the imaginary part of the self-energy. In 
this article we offer a real time critique of this relation by studying the real time evolution 
and relaxation of linearized perturbations, and point out the following observations: 
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1. The damping rate is identical to the imaginary part of the self energy only when the 
inflaton (the particle interpolated by the order parameter) is a resonance with T being 
its width. In the complex frequency plane this corresponds to a pole in the second 
(unphysical) Riemann sheet. Even in this case, this damping rate corresponds to 
the relaxation of the expectation value of the scalar field and is exponential for some 
time regime. However, eventually relaxation continues with a power law tail. When 
the imaginary part is zero on-shell and the one particle pole is below the multiparticle 
thresholds, relaxation is given by a power law and no “damping rate” can be associated 
with the imaginary part of the self energy. Moreover, the relation T ps ImE (ou,k)/ou 
only holds in the linear approximation around equilibrium. 


2. Even when relaxation of the expectation value of the scalar field is exponential, this 
damping rate determines the approach to equilibrium of the expectation value of the 
scalar field but it cannot be immediately inferred that the same time scale describes 
thermalization, i.e. the approach to a Bose Einstein distribution of an initial off- 
equilibrium distribution. Thermalization is a rather different process and has to be 
described with a Boltzmann equation with a collision term. Thus, rather than inter¬ 
preting this time scale as a thermalization scale, we interpret it as a relaxation scale 
for the expectation value. 


3. We consider the non-linear quantum field evolution (non-linear relaxation) incorpo¬ 
rating the order parameter into the effective mass in a self-consistent way. [The self- 
consistency requirement makes the evolution equations non-linear]. We study the 
inflaton evolution in real time coupled both to light scalars and fermions, providing 
exhaustive numerical results. We find that the non-linear relaxation time scales can 
be much shorter than those predicted by linear relaxation. In addition, the particles 
are produced with a momentum distribution that is very far from thermal and skewed 
towards low momentum. In the case of fermions we provide numerical evidence for 
the phenomenon of Pauli blocking that hinders dissipation and production of fcrmion- 
antifcrmion pairs. This phenomenon is very similar to that found by Kluger and 
collaborators 


23fl in their studies of fermionic back-reaction in the presence of strong 


electric fields. 


Both time scales will have to be understood in detail to provide a quantitative and reliable 
estimate of the reheating temperature. 

Although we are ultimately interested in describing reheating and thermalization in 
Friedmann-Robertson-Walker cosmologies, in this article we will work in Minkowski space 
assuming that the relaxation time scales are much shorter than the expansion time scale. 
Furthermore, we want to isolate dissipative effects arising from particle production from 
those resulting from the red-shift in an expanding cosmology. 

A short account of the results presented in this paper was reported in ref. 
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In section II, we introduce the general model that we propose to study and summarize 
the necessary ingredients of non-equilibrium held theory to provide the reader with the 
essentials needed to reproduce our calculations, as these do not seem to be part of the 
standard techniques. In section III, we study linear relaxation in real time and elucidate 
the role of the imaginary parts of self energies, their interpretation in real time and criticize 
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the identification with “thermalization rates”. In section IV we study the process of non¬ 
linear relaxation, and particle production both for scalars and fermions, providing exhaustive 
numerical results. Finally we conclude in section V with a discussion of the implications of 
our results and provide a guide to obtain estimates for the reheating temperature in model 
theories. An Appendix is devoted to a pedagogical exercise in non-equilibrium field theory. 


II. THE MODEL AND THE METHODS 

We consider the simplest model |2].|d] where the inflaton field $ couples to a scalar a and 
to a fermion field 0. That is, 

£ = (S 2 + m| + go 2 ) 4 - 4 4 - iff (9 2 + m 2 ) o - A a ' + - yl'jy . 

( 2 . 1 ) 

The case m a , will be of particular relevance, since in this case there are open 

decay channels for the inflaton. 

We will investigate the scalar and fermionic couplings independently. Although the 
situation for reheating corresponds to (almost) zero temperature, we will study the case of 
linear relaxation at finite temperature. The reason for this is that finite temperature effects 
reflect the Bose enhancement and Pauli blocking factors that appear whenever there are 
excitations in the medium. These contributions from the medium will allow us to identify 
similar physical features in the case of non-linear relaxation. 


A. Linear Relaxation: Amplitude and Perturbative expansion 


To explore the behavior of the inflaton within the linear regime, we first use the tadpole 
method to obtain the equation of motion (|A3|) (see the Appendix for details on how this 
method is actually applied and reference |25j for an alternative implementation). The next 
step is to linearize this equation in the inflaton zero mode amplitude and use this to study 
the relaxational dynamics of the inflaton. We should note that this amplitude expansion is a 
priori different from the standard perturbative expansion in the relevant coupling constant. 
Later in this work, we will compare the results obtained here with results obtained through 
a self-consistent, non-perturbative resummation both in the coupling constants and the held 
amplitude. 

We arrange the initial conditions to be such that the fields <f>, a, 0 start to interact at a 
time that we choose t — 0. This can be accomplished by making the coupling “constants” 
time dependent, i.e. zero for t < 0 and different from zero for t > 0. To perform the 
calculations, we will need the non-equilibrium Green’s functions and Feynman rules. Since 
the non-equilibrium generating functionals involve a forward and backward time contour 
|8].[22|,|26|] , the number of vertices is doubled. Those in which all the fields are on the forward 


branch (fields labeled by (+)) are the usual interaction vertices, while those in which the 
fields are on the backward branch (fields labeled by (—)) have the opposite sign. The 
combinatoric factors are the same as in usual field theory. The spatial Fourier transform of 
the necessary finite (initial) temperature propagators are: 
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• Bosonic Propagators 


Gi + (t, t') = G>(t, t')Q(t - f) + G<(t, - t ) , 

Gr(t,t ') = G>(t,t')0(t’ -t) + Gk(t,t')Q(t - t') , 

Gi-(t,t') = -G<(t,t'), 

G^(t,t>) = -G>(t,t'), (2.2) 

G%(t,lf) =i I d 3 x e~ it3 <$(£,t)$(0,t')> 

= ^ {[l + , 

G^(M') =* J d 3 xe~ itx {$(0,t')$(x,t)) 

= ^ {[! + n & (a; fc )]e^- 0 + n^e^- 0 } , 
o; fc = \[W^rn 2 , n 6 (o; fc ) = —- , 


where m is the mass of the boson. 

• Fermionic Propagators (Zero chemical potential) 

Si+(M') = Sf(MW - 0 + - t ) , 

srM = s2(t,t)ev -t) + ss(t,t')o(t - 1 '), 

Sj-(t,t') = — i J d 3 x e~^' x (^(x,t)i}(0,t')) ( 2 . 3 ) 

= [e _,Wfc(t_t,) (^ + m v ,)( 1 - 71/(0;*,)) + e* aJfc(i_i ' ) 7 0 (/: - m^)'y Q n f (a;*)] , 

S~(t,t') =ij d 3 x e~ l ^' x (- 0 ( 0 , t')il>(x, t)) 

= + e -^(*-*') 7o (^ - m V) ) 7o (l - n/(a; fc ))l , 

zuk L J 

+ = e/? ^ + x • 


In the linear amplitude approximation, corresponding to linear relaxation, we find in all 
cases the following form of the equation of motion for the expectation value (see Appendix) 

poo 

Spit) + nl 5 P it) + / Kp(t - 0 %(f') dt' = 0 , (2.4) 

J 0 

A>(t - t') = E r> p(t - f') 0 ^ - t') , 

where we have imposed as boundary conditions that the inflaton and the other fields are 
coupled at time t = 0 but uncoupled for previous times, and introduced 
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(2.5) 


8p(t) = j d 3 x e l V' x 4>(x, t ), til = p 2 + m% + 5m{T) , 


where Sm(T) is the time independent (but temperature dependent) contribution from tad¬ 
pole diagrams. These contributions renormalize the mass and introduce a temperature 
dependent effective mass and will be specified later in each particular case. The quantity 
E r ,p(t — t')0(t — t') is the retarded self-energy. It will be computed to dominant order in the 
couplings for both fermions and bosons. Although the self-energy at finite temperature has 
been computed before in the literature P7| , |28| , we differ from previous treatments in that 
we perform our calculations directly in real time, this allows us to study real time relaxation 
as an initial condition problem. 

Eq. ([2.4|) can be solved by Laplace transform. Define 

/*oo 

<pp(s) = / e~ st 8p(t) dt . (2.6) 

J 0 


Then, eq. (|2.4|) becomes 


s 2 w(s) - s<5jr(0) - 4(0) + SI} <Pp(a) + ipp(a) Ep(s) = 0 , 


(2.7) 


with Sp(s) the Laplace transform of E r< p(t). 

For computational purposes, it can be shown that at zero temperature each graph of E p(s) 
is exactly equal to the corresponding graph of the zero-temperature equilibrium Euclidean 
quantum held theory, s being the time component of the Euclidean four momentum. 

In general E p(s) can be written as a dispersion integral in terms of the spectral density 

p(Po,P,T ) 



2 PoP(Po,P,T) 

s 2 +pl 


dp 0 ■ 


The imaginary part of the self-energy is found to be 


( 2 . 8 ) 


Im E p(s = iu ± 0 + ) = ±Y,jp(u) , 

Ejp’(a’) = vr sign (cu) \p(\u\,p, T) — p(—\u\,p, T)\ . (2.9) 


The presence of sign (cu) in the above expression characterizes the 
Let us choose 5^(0) = 8i , 5^(0) = 0 for simplicity. We get from 


retarded self-energy, 
eq. 0 


P p( ) 8i 


s 

s 2 + tip + E p(s) 


( 2 . 10 ) 


The Laplace transform can be inverted through the formula 


r+ 

8p(t) = / 

J —i 


+ZOO+6 (j a 

e st pp(s 

—200+6 


27t i ’ 


( 2 . 11 ) 


where e is a positive real constant (Bromwich countour). Thus we need to understand the 
singularities of [s 2 + tii+ Ep{s)]+ We now consider the following cases: 


• The inhaton potential admits spontaneous symmetry breaking (SSB), and is only 
coupled to lighter scalar fields. 
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• The inflaton is coupled to fermions only. 

We will also consider the subcases where the initial temperature is taken to be zero, as 
would be appropriate in the case of evolution in the post inflationary universe, as well as the 
situation where the initial temperature is non-zero, which would be relevant to the situation 
of a scalar held starting in an initial (non-equilibrium) but thermal state and evolving out 
of it. 


B. SSB with Coupling to Lighter Scalars Only 


If m| = —p 2 < 0, then the new minimum is at <f>o = \j 6/i 2 /A$, and we write <f )± (x, t) = 
$0 + t ) + y ± (^, t). The masses are now shifted to 


M 2 = 2/i 2 , 


Ml = ml + g 2 & 


( 2 . 12 ) 


The contribution from the quartic inflaton self-coupling has been studied previously III 
thus it will not be repeated here. 

The tadpole correction to the mass in eq. ( [2.5|) is given by 


8M(T) = -ig j — g J 


d 3 k 1 + 2 n b {u k ) 
M 3 2o4 


u k = V k 2 + ill 2 


This mass renormalization does not influence the dynamics. The retarded self energy is 
found to be at order g 2 , 


Kp(t - t') = 2 ig 2 $l 


d 3 k 

(2^ 


++ - + - 0 - Gf/t - t')Gi +p> (t - 0] ■ (2.13) 


With the non-equilibrium Green’s functions defined above, we find 

d 3 k 1 


KAt - f> ) = - / 


(2tt) 3 2 

-2 n b (uj k ) sin[(o;^ - u n+jf ){t - 1)]} 


{[1 + 2 n b (uj k )] sin[(o;^ + (t - if)] 


(2.14) 


The Laplace transform can be written as a dispersion integral in terms of the bosonic spectral 
density p b (p 0 ,p,T) (see eq. (|2.8| )). We have to one loop level: 


Pb{Po,P,T ) = 2 g $ 0 


d 3 k 


d 3 k' 


{2n) 3 5 3 {p— k — k') 


(2vr) 3 2 uj k J (2vr) 3 2 u k ' 
x [5(p 0 — uj k — o;fe/)(l + 2n b (uj k )) - 8(p 0 - uj k +u fc -)2n t (+)] • 


The imaginary part of the self-energy is given by eq. (ED- 

We will analyze only the case of a spatially constant order parameter corresponding to 
p — 0 because we will later compare with the case of non-linear relaxation which we only 
study for a homogeneous (translational invariant) expectation value. In this case the spectral 
density can be written as 
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Pb{p 0 , 0, T) = 


1 + 2 n b (j 


p t (p o ,0,0) 


(2.15) 


The spectral density pb(p 0 , 0, 0) is a Lorentz scalar and is proportional to the decay rate of 
the boson $ into two cr particles. It is a straightforward exercise to find 


Pb(Po, 0,0) = 


g 2 H 

87r 2 


4 Ml 
pi 


e (Pc - 4Mj) . 


(2.16) 


Clearly Eg(s) has a logarithmic divergence, independent of s and T. We choose to subtract 
this divergence at s = 0 and absorb the subtraction into a further temperature dependent 
renormalization of M. In order not to clutter notation, from now on we will refer to M as 
the fully renormalized mass including the above subtraction, and to E(s) = E^(s) — Eg(0). 
From equations ( \2 . 8| . f2 . 1 5|. |2.16| ) we hnd that the self-energy has an imaginary part above 


the two cr-particles threshold given by eq. (|2.9|) with 


E/M = 


sW 

8n 


1 - 


4 Ml 


c 0- 


1 + 2 iib 


0 (a; 2 — 4M 2 ^) sign (a;) 


(2.17) 


This expression is recognized as the imaginary part of the retarded self-energy (determined 
by the sign (a;)) and shows the usual Bose enhancement factor |29| . 


1 . Zero Temperature 


In this case from |2.14| we hnd 


S (s) = 


47r 2 


4 M 2 

1 + —ArgTh 


1 _l 4M| 
1 + s 2 


(2.18) 


In order to compute the inverse Laplace transform through eq. (|2.11|) we must first study 
the analytic structure of <p(s) in the s-plane. <p(s) has poles at the zeroes of 


s 2 + M 2 + E(s) = 0 


(2.19) 


These correspond to 4)-one-particle states with the mass including one-loop radiative cor¬ 
rections. At zeroth order the poles are purely imaginary 

s± = ±iM . (2.20) 


To hnd the one-loop correction, we set 

( 2 . 21 ) 

( 2 . 22 ) 

That is, 


s_|_ = %M + r 

and similarly for s_. Inserting this in eq. ( |2.19| ) yields to order g 2 

2iMr + E (iM) = 0 . 
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(2.23) 


r = 


.E(iAf) 

1 2 M 


When M < 2 M a , E (iM) is real (see eq. (|2.18|) ) and eq. ( p.23| ) gives a real correction to the 
<h mass 


S-i- = = ± i 


M 


9 2 ®t 


8 t r 2 M \ V M 2 


4M 2 4 1 

1 ArgTh ■ 


AM 2 

M 2 


(2.24) 


The Laplace transform <p(s) also exibits a cut along the imaginary axis starting at s = iu = 
±2 iM a . For s in the first Riemann sheet (physical sheet) we obtain 


^physical {iu i 0 ) — E/j(cn) A iE/(cn) , 
with E r and E/ both real and given by 


u > 2 M a , 


(2.25) 


Sr(cu) — 


47T 2 


' 1 - 


4M 2 




ArgTh \ 1 — 


4M 2 


S r(u) = 


l 


g' 2 ®l 

8tt 


4 Ml 

'1-^ > 0 . 




(2.26) 


We can now proceed to compute the inverse Laplace transform ( |2.11| ) by deforming the 
contour. 


m 


St cos M 0 t | 2 Si f°° uT,i(u) cosut du 

+ V hu a [a7 2 -M 2 -EHM] 2 + E 7 H 2 ‘ 


(2.27) 


For large time M a t 1 the integral over the cut is dominated by the endpoint u = 2 M a 
and goes to zero as 


S cn t{t) ~ 


Sj^Mlg^l cos(2 M a t+?f) 

4tt 2 (M 2 - 4M 2 + ^) 2 (M ff f)l 


(2.28) 


The f -3 / 2 is completely determined by the behavior of the spectral density at threshold. 

The situation changes drastically for M > 2 M a . In such case the T-particle is unstable 
and thus E (iM) becomes complex and its value depends from which side of the cut we 
approach the imaginary axis. Now the solution for the pole will be complex and r will 
acquire a real part. Due to the discontinuity in E/ across the two-particle cut, eq. (j2.23| ) 
must be written carefully as 


r = 


2 M 


E (iM + Re [r]) , 


(2.29) 


where the (small) real correction inside the argument will determine on which side of the 
cut the solution resides. Then the real part of r should satisfy 


Rer = 


2 M 


Im E (iM + 0 + sign (Re [r])) . 


(2.30) 
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From eq. (|2.17|) we see that ImE p h ys icai is negative for sign (Re [r]) < 0 and positive for 
sign (Re [r]) > 0. Therefore eq. ( [2.30D has no solution in the physical sheet. 

The analytic continuation of £ into the second Riemann sheet is such that 


3C 


y. n ( 


-u rr 




and we hnd the solution 

1 


Rer = 


2 M 


MM) = 


Im r = 


2 M 


MM) = 



< 0 


ArgTh \ 1 — 


4 Ml 


For M M„ 


Im r = 


g 2 ® o 


M 


8vr 2 M V l0g M„ 1 


(2.32) 


| Re [r]| coincides with the decay rate 4> —> 2cr (as it must be) which is the rate per unit time 
to produce a particles. The poles s± move off into the second sheet when M becomes larger 
than 2 M a as expected . 


Thus when we compute the inverse Laplace transform (|2.11|) in the unstable case we are 
left with the integral over the cuts since both poles are in the second Riemann sheet and we 
hnd 


m 


2 Si r°° ojTji(uj) cos oot du 

i r J2M a [a; 2 — M 2 — Sr(o;)] 2 + £/(u;) 2 


(2.33) 


Since now M is inside the integration region, for weak coupling there is a narrow resonance 
at uj ~ M. Thus for weak coupling it takes Breit-Wigner form and we hnd to a very good 
approximation, 


5(t) ~ Si A e rt1/2 cos (Mt + a) , T -C M , 


where 


dMM) g 2 ® 2 0 I 4 Ml = dMM) 

dM 2 ’ 8ttmV M 2 ’ ° dM 2 


For M M a , 


A 


1 + 


g 2 ®l 

8vr 2 M 2 


+ o 



g 2 ^lMl 

4ttM 4 


+ o 



(2.34) 


(2.35) 


(2.36) 


The Breit-Wigner approximation, however, is valid only for times < ln(T /M a )\ for longer 
times the fall oh is with a power law determined by the spectral density at threshold 
as before. 
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2. Non-Zero Temperature 


The physical mass gets a finite temperature dependent correction from the tadpole con¬ 
tribution 5M(T) given by 


[°° dk k 2 1 

5M(T) =g —- =- , 

J0 2tT 2 Up + M 2 /3y/k 2 +M% 


(2.37) 


For small (3 (T M a ) this correction takes the form 

T M„T Ml 


w < T > = »{S-^ + # l 0 Uy> + o < rC, >}- 


(2.38) 


We will assume that gT 2 -C M 2 so that we are in the perturbative regime, since otherwise 
hard thermal loops must be resummed [3TIJ32], a task beyond the scope of this article. 

The imaginary part of the self-energy is given in eq. (|2.17|) and the real part can be 
obtained from the dispersion integral ( [2.8p using ( |2.15| , [2.1(j| ) . The real part of the self-energy 
is difficult to compute for arbitrary temperature, and we just quote its large T behavior: 


M T 

S (s,P) = g*** a 


° ITS 2 



s 2 

1- 

1 

+ 

s; 

q to 

1_ 


+ O(T 0 ) 


(2.39) 


It is interesting to compute E (t,/3) in configuration space for large T. We have by inverse 
Laplace transform 


/ +ioo+e We 

e s *E( S ,/3) —. 

-ioo+e 2m 

Upon contour deformation we find 


n 2 ^>xT r°° dr ,_ 

E (t,(3) =--— / — Vx 2 — 1 sin(2 M a xt) . 

m J i x 


(2.40) 


(2.41) 


This function is obviously not concentrated at t — 0. It can be related to a J\ Bessel 
function as 


S (t,/3) = 


2vr 


am r 1 

1-- t — 2 M„ / dx 


TT 


t 


X 


1) Ji(2 M a x) 


(2.42) 


We hnd for small and for large t 

EM) f =° - 


g 2 <b 2 T 


2tt 

a\ g^lT sin(2 M a t - tt/ 4) h , 

Ht ’ p) ~ ~2 v (j^t ym + > 


i_LL i + 0 (i 2 ) 

7T 


The behaviour of the kernel E(t, (3 ) shows that it cannot be approximated by a phenomeno¬ 
logical term T ^ even for high T. (As shown in ref. |L| this is not the case either for T = 0). 
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We can now repeat the analysis of the previous section for both cases M < 2 M a and 
M > 2 M a . In the first case, the one particle pole is below the two cr-particles threshold with 
a (small) finite temperature correction since we are restricted to the perturbative regime in 
which gT 2 -C M. The long time behavior will be oscillatory with the frequency correspond¬ 
ing to the one particle pole plus long-time power law tails similar to the zero temperature 
case. The second case is more interesting, since now the scalar “inflaton” is unstable, and 
the pole moves off into the second Riemann sheet. In the physical sheet there is a resonance 
with a finite temperature width given by 


T(T) = T(0) 


1 + 2 Tib 



T(0) 


I 4 Ml 

8tiM V M 2 


(2.43) 


The Bose enhancement factor increases the rate and therefore enhances dissipation via the 
production of particles. Although this factor arises from the thermal distribution, we expect 
in general that whenever there are bosonic excitations present, the relaxation rate will be 
enhanced as a consequence of bose statistics, independently of whether these excitations are 
thermally distributed. 


C. Unbroken symmetry with coupling to light 


scaiars 


In the unbroken symmetry case is above all thresholds and hence $ is a stable 
particle. The order parameter is again given by a formula like eq. (|2.27|) except that now 
the integration starts at u = M + 2 M a . 

The first perturbative contribution to the kernel E (t — t') in the inflaton equation of 
motion (|2.4|) are now the two loop order graphs usually called “setting sun”. Since these 
graphs are quite complicated, we only perform the two loop computation at zero temper¬ 
ature, for which we can exploit the relationship with usual Euclidean field theory. Even 
with this simplification the computation is hard, and we limit ourselves to the evaluation of 
the imaginary part of the retarded self-energy near the branch point. We have seen below 
eq. (|2.27|) how this actually determines the long time behaviour of the order parameter. In 
general, if E/(u; —> a; t h re shoid) vanishes as (a; — threshold) 0 , then 5 cu t(f) decays as U 1_a for 
large times. 

We find for the two loops self-energy 


E/(u; —> M + 2 M a ) ~ 


2g 2 7T 2 M a \fM r 2 


(4vr) 4 (M + 2M a ) 7 / 2 


[u 2 - (M + 2 M a yy • 


(2.44) 


This yields for 5 C ut(^) a power decay t 3 for long times. 
For massless a we find 


E/(u; —> M) ~ 


g 2 n 


■\co 2 -M 2 ) 3 , 


3M 4 (4tt) 


(2.45) 


yielding a 5 C ut(^) decaying as t -4 for long times. In summary, in the unbroken symmetry 
case 5{t) is given by the one-particle term oscillating with frequency M 0 plus the power 
damped cut contribution 5 C ut(£)- 
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D. Inflaton coupled to Fermions 


Fermions can be treated similarly in the broken and unbroken symmetry phase. To treat 
both cases on equal footing we now define M as mass of the inflaton (scalar) held and m as 
the fermion mass in either case. Using the Feynman rules described in the first section, we 
obtain to one-loop level the kernel 


Kp(t - t ') = iy 2 
K,p(t - t') = iy 2 


d 3 k 

(2^)3 

d 3 k 
(2vr) 3 


Tr 

Tr 


t')iSi+ f (t', t ) - iS£(t, t')iS~_p(t', t ) 
iS? (t, t ) - , t ) 


We now concentrate on the homogeneous case p = 0. Using the fermionic Green’s functions 
given in the first section it is straightforward to find the fermionic spectral density to be 
used in eq.([2.8D 


Pf(Po,T) = ^pl 


1-2 n 


f 


1 - 


4m 2 

pI 


0 [p 2 0 — 4m 2 ) 


From eq.(2.9), the imaginary part of the self energy is then 


v / x _ v 2 2 
E,(w) “ 8?" 


1 - 


4 m 2 


a r 


1 - 2 n f ( “ 


0 ^u; 2 — 4 m‘ 


The hnite temperature factor rehects the Pauli blocking term |^9[. It is clear that the zero 
temperature part of £(s) diverges quadratically and that two subtractions are needed. The 
first one is independent of s and contributes a (quadratically divergent) mass renormaliza¬ 
tion. The second one is logarithmic divergent and consists of an s independent term that 
adds to the mass renormalization and another proportional to s 2 that will be absorbed in 
wave function renormalization. We choose to subtract at zero temperature and at an arbi¬ 
trary scale k. The Laplace transform for the zero momentum component of the equation of 


motion (|2.10|) becomes 


Vis) = 


Si s 


s2 S + S ln h + VisCD + v 2 n (s, t, K ) ’ 


(2.46) 


where M 2 R contains the mass renormalization and II is the twice subtracted kernel, which 
depends on the renormalization scale. Defining 


V R (k) = zjy , 

M r (T, k) — Zj M iR (T) , 
(Pr{s, k ) = Z^ 1 <p(s) , 


we finally obtain the renormalized Laplace transform 
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Si s 


s 2 + M 2 R (T,K) + y 2 R (K) II (s,T,k) 


(2.47) 


The inverse Laplace transform is obtained as in equation (|2.11|) . The result will be the 
function 5 R $(t, k) which is not a renormalization group invariant. It is clear, however, from 
the renormalization prescriptions described above that ratios of the amplitude at different 
times, such as S R $(t, k)/Sro(0, k) are renormalization group invariant. Now the analysis 
can proceed as in the previous section. To obtain the inverse Laplace transform we must 
recognize the singularities in ( |2.47p . If Mr < 2m ( m is the fermion mass in the loop), there 
is a one particle pole (with strength different from one because we decided to renormalize 
off-shell) and a two-fermion cut at s 2 = —4m 2 . At long times the amplitude will oscillate 
with an oscillation period given by the position of the pole, which is perturbatively close to 
M r . The contribution of the cut falls-off at long times as f~ 5 / 2 and is determined by the 
behavior of the spectral density near the two-fermion threshold. 

More interesting is the case in which Mr > 2m. As in the bosonic case, the pole moves 
off the physical sheet into the second sheet. In the physical sheet the spectral density at 
weak coupling features a sharp peak at Mr + 0(y 2 ) with width 


T(T) = T(0) 


1 - 2 n 


M 


R 


r <°> = k M « 


1 - 


4 m 2 


(2.48) 


A Breit-Wigner approximation predicts exponential relaxation but eventually at long times 
(an estimate similar to the bosonic case) a power law relaxation t~ 5 / 2 ensues, completely 
determined by the spectral density at threshold. 

Pauli blocking makes the resonance narrower and the lifetime of the decaying particle 
longer. The interpretation of this phenomenon is simple. In the thermal bath, fermionic 
excited states are filled with the Fermi-Dirac distribution. In order for the scalar field 
to decay, it must create a fermion-antifermion pair, ffowever, at finite temperature, the 
available states are already filled with thermal excitations and because of the Pauli exclusion 
principle, are not available. At infinite temperature (and zero chemical potential), each 
fermion and antifermion state are populated with occupation 1/2 per spin degree of freedom; 
in this limit the decay rate goes to zero, the lifetime to infinity and the bosonic particle simply 
cannot decay because there are no states available to decay into. Even at zero temperature, 
but in a situation in which excited states are occupied, dissipative processes mediated by the 
production of fermion-antifermion pairs will be hindered by the Pauli exclusion principle, 
since states will already be occupied and no longer available in the particle production 
process. In highly excited states, we expect damping via production of fermion pairs to be 
strongly suppressed by Pauli blocking. This phenomenon has been seen numerically in the 
case of fermion pair production in strong electric fields |^3| and will be seen numerically in 
the non-linear relaxation case later. 


E. Thermalization or Relaxation? 

From the analysis presented above, the following conclusions for the linear regime become 
very clear: 
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1. The imaginary part of the self-energy only determines a relaxation rate in the case of 
a resonance, that is when there is an imaginary part on-shell for the external particle 
and the (quasi) particle pole moves off the physical sheet into the second (unphysical) 
sheet. In this case a Breit-Wigner approximation describes the exponential relaxation 
for a long time, but eventually the amplitude falls-off with a power law in time, with 
the power determined by the behavior of the spectral density at threshold. When the 
on-shell pole is below the two-particle threshold, the imaginary part does not translate 
to a damping rate. Relaxation is described by a power law with an asymptotic behavior 
completely determined by the position and residue of the pole. 

2. Even in the case of a resonance and exponential damping, the “damping rate” T 
describes exponential relaxation for the expectation value of the scalar field. The issue 
of thermalization is completely different. Thermalization corresponds to the time 
evolution of the (quasi) particle distribution function towards a thermal distribution 
which, in principle, has nothing to do with the relaxation of the expectation value of 
the field. 

An alternative way to look at thermalization is as a process of momentum and energy 
transfer. Thus, the thermalization rate should be identified with the rate of energy and 
momentum transfer which is not necessarily related to the relaxation rate of the expectation 
value of the field. In order to understand thermalization, a collisional (quantum) Boltzmann 
equation must be set up. Although in the Born approximation the collision term includes 
the scattering cross section that is related to the decay rate, the solution to the Boltzmann 
equation implies a resummation that is quite different from the resummation of the Dyson 
series for the propagator. If the initial distribution is very far from thermal and many 
collisions are necessary for thermalization, the relaxation and thermalization time scales 
may be widely different and in principle unrelated. Thus we insist that the “damping rate” 
obtained from the imaginary part of the self-energy on shell must be interpreted as the 
relaxation rate for the expectation value of the scalar field and in principle not with the 
thermalization rate of the particle distribution function. Moreover, this relaxation rate 
holds in the linear regime (small field amplitude). 


III. NON-LINEAR RELAXATION 


In this section we study the equation of motion for a homogeneous order parameter 
beyond the linear approximation. This is achieved as follows [|T3|. We write the scalar 
inflaton field as 


<h ± (x, t) = (p(t) + x ± (x, t ) , (3.1) 

with «p(t) the expectation value in the non-equilibrium density matrix and < x ± (x, t) >= 0, 
we consider that < a ± >= 0 (this is a consistent assumption). The non-equilibrium path 
integral requires the Lagrangian density 


£0;x + w + ;x 
£0;x + w + ) = 


W ) = >C(0;x + ;ct + ) - £(0;x ;a ), 

1 r 


i (d,X + ) 2 - M‘\t) (x + ) 2 


2 L 


X 


A 


<t> + ml 0 + - 0 3 

o 


(x + ) 3 

6 


A 

4! 


(x + ) 4 
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^ [(<V + ) 2 -m 2 (t) (cx + ) 2 ] - gcpx + (° + ? ~ | (X + ) 2 {(J + f + - yx + )^ + , 

r2 ' -' '> ^ 


M 2 (t) = m| +-</> 2 (t) , m 2 (t) = m 2 + g </> 2 (t) , m^(£) = + 7/0(t) 


The difference with the linear relaxation case is that we now incorporate 4>(t) in the definition 
of the time dependent masses. The equations of motion are obtained as in the previous 
sections, by treating the linear, cubic and quartic terms as perturbations. The necessary 
Green’s functions are constructed from the homogeneous solutions of the quadratic forms. 
We again treat the cases where the inflaton is coupled to scalars or fermions separately. 


A. Inflaton Coupled to Scalars Only 


The Green’s functions for the scalars are obtained from the mode equations that solve 


+ k 2 + M 2 (t) U k (t) — 0 , 


£4(0) = 1 , £4(0) = -iW k = -i\jW + M 2 (0) , 

v k (t) = 0 , 


^ + P + m?{t) 


Vfc(O) = 1 , 14(0) = -iw k = —i\]k 2 + m 2 (0) , 


(3.2) 


(3,3) 


The initial conditions on the mode functions (T2, T3) correspond to positive frequency 
solutions at the initial time. In terms of these mode functions, the Green’s functions are 

01 


G? lfc (M0 = [(1 +n b (W k )) U k (t)U* k (t') + n b (W k ) U* k (t)U k (t')] , 


G> k (t,t') = — [(1 +n b {w k )) V k (t)V k (t') + n b (w k ) V k [t)V k {t')} , 

G< k (t,t') = G> k (t',t) 

and the rest of the Green’s functions are given by the relations in equations (|2.2| ). The 
Green’s functions above correspond to the situation in which the initial density matrix is in 
equilibrium for the positive and negative frequency modes at time t = 0. This is determined 
by the initial conditions on the mode functions above at an initial temperature T = 1 //3. 
However, for the rest of the analysis we will take T = 0. We see that 


< X 2 (x,t) >= J 

< cr 2 (x, t ) >= J 


d 3 k \U k (t )\ 2 
(2vr) 3 2H4 
d 3 k \V k (t)\ 2 
(27r) 3 2 w k 


coth 

coth 


' pw k 

2 

fdw k 


Finally the equation of motion to one-loop order for the expectation value is 
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(3.4) 


4>(t) + m| j + 2 4>(t) < X 2 (x, t) > +9 < a-(x, t)>= 0 

o 2 


If one wants to solve the equation of motion (|3.4|) to order ( h ) one would expand 0 in a 
power series in h and only keep the zeroth order term in the mode equations (|3.2| , |3.3| ). As 
was observed previously [[H|], such an expansion will result in secular terms and becomes 
unreliable at long times. A resummation is necessary to capture the long time behavior 
consistently. We will perform a non-perturbative resummation of the one-loop terms by 
solving the set of equations (fO - |3.4|) with the full value of 0 in the mode equations. 


This then becomes a set of coupled non-linear integro-differential equations that provide 
a non-perturbative resummation of select one-loop terms as can be seen by looking at a 
perturbative expansion of the solution to these equations. 

We want to emphasize this point. By incorporating the full value of 0 in the evolu¬ 
tion equation for the mode functions we are incorporating back-reaction effects. If only the 
classical evolution of 0 is used in the mode equations, we would have parametric resonant 
amplification since the effective frequencies for the mode functions are periodic due to the 
fact that the classical solution is periodic with constant amplitude. This leads to particle 
production that never shuts-off. However particle production leads to damping in the evo¬ 
lution of 0. Introducing this damped evolution in the mode equations leads to a behavior 
rather different from parametric resonance: as the evolution of 0 is damped, the amplitude 
becomes smaller and particle production should diminish and eventually stop. This was 
found to be the situation in the case of the scalar held with self-interaction |I3| . 

Clearly the integration of this set of equations will have to be done numerically. Because 
the self-interacting scalar case was already studied before B we will now concentrate on 
the interaction with the scalar holds in which thresholds are present. But before we do this 
we must understand the renormalization aspects of this system of equations. 

The large-fc behavior of the mode functions can be understood from a WKB analysis, 
the details of which had already been discussed in [||,|T3j. We find the following divergence 
structure doing this 


d 3 k \U k (t )\ 2 _ A 1 

(2vr) 3 2 W k ~ 8^ “ 8^ 
r A d 3 k \V k (t)\ 2 A _ 1 

(27t) 3 2 w k 87 r 2 87 t 2 


m 


+ 'y0 2 ( i O j l n — + Fi{t, k) , 


(' m l,b + 9 b(p 2 (t )) In ^ + F 2 (t, k) , 


where A is an upper momentum cut-off and n and arbitrary renormalization scale. The 
subscript b refers to bare quantities and the quantities F 1)2 (C ft) are finite in the limit A —> 00 . 
We can now read the mass and coupling constant renormalizations 


m i >,r ~ 


A_r — Xb — 



8n 2 


2 1 2^1^- 

Y m $,b + 9b m a b 1 In — , 


We introduce a further, finite renormalization, by subtracting the functions A lj2 (t = 0, hi), 
absorbing this subtraction into a (finite) renormalization of the mass. After these renor¬ 
malizations, we finally arrive at the renormalized set of evolution equations in terms of 
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renormalized quantities. We drop the subscript R for renormalized quantities to avoid clut¬ 
tering the notation, but with the understanding that all quantities are renormalized at the 
scale k: 


4>{t) + rn% <f>(t) + ^ </> 3 (t) + ^ <j>(t) 


f A d 3 k \U k (t)\ 2 -l 


(2tt) 


2 W k 


+ 9 4>{t) 


'A 2 


87T 2 


+ 9 2 )m 0) 


A 

In — = 0 . 

K 


f A d 3 k \V k (t)\ 2 -l 

J (27r) 3 2w k 


(3.5) 


To this order, we can replace masses and couplings by their renormalized values in the 
equations for the mode functions (|3.2| , |3.3|) . We will now chose the renormalization scale 
to be k = rn<i, so as to have only one scale in the problem, which makes the numerical 
evaluation easier. 

The renormalized equation of motion ( |3.5|) may be written without reference to the 
cutoff A which in the end must be taken to infinity, however, numerically the /c-integrals 
must be calculated with an upper momentum cutoff. One must ensure that this numerical 
cutoff, to be identified with A in the evolution equation, be much larger than the masses and 
amplitudes of the field for the integrals to reach their asymptotics and the cutoff dependence 
to dissapear. We have been careful to make exhaustive checks that the final results were 
insensitive (to the working accuracy) to the cutoff which was typically chosen to be A & 
100 This implies that the order of magnitude of the error is about (^) 2 — 10~ 4 . 

A provision must be made for the initial conditions on the mode functions solution 
of ( |3.2| . |3.3| ) in the broken symmetry case. In this case the tree level squared mass is 
negative m\ = — /i 2 and for initial values of the expectation value 0 2 (O) < 2/i 2 /A the initial 
configuration is below the classical “spinodal” and imaginary frequencies lead to the spinodal 
instabilities. We are not interested in this article in studying the time evolution of these 
instabilities but on the issues of non-linear relaxation. There are two ways to avoid the 
complex frequencies associated with these instabilities in the initial conditions: (i) one can 
choose an initial value </> 2 (0) > 2/j 2 /X or (ii) one can impose that the initial frequencies 
correspond to a positive mass term. This choice corresponds to preparing an initial gaussian 
density matrix of the modes with this given mass, under time evolution this packet spreads 
in function space and the time dependence of the width reflects the (quantum and thermal) 
fluctuations. In our numerical analysis we chose the later possibility with the positive mass 
squared given by the absolute value of the (negative) mass squared in the Lagrangian. 

We introduce now the following dimensionless variables: 


V(r) 


\L\ 2 1^)’ 

V 6 \m%\ 

9 = 

W q 

~\ 

2 |A/2(0) 

l"4l 

W q 


k 

\m*\ 


r = |m$| t , 



m 2 ( 0) 

\m%\ 


In terms of which the evolution equation (|3.5|) becomes 


v( T ) + v(t) + V 3 ( t ) + ^( r ) 


A 

n<j> | 2 


q dq 


\U q {r)\‘ 


^ il + %) " (T) “ ’' 2(0) 


W n 


In 


A 




1 + 

= 0, 


'^ q 2 dq iwi: 

Wn 


(3.6) 
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where we chose the renormalization scale k = 


1. Particle Production 


As in any time dependent situation, the concept of particle is ambiguous and has to be 
specified with respect to some particular state. We choose that state as the equilibrium 
ensemble at the initial time t — 0. The initial condition on the mode functions for the 
fluctuations Q3.2| . |3.3|) naturally determine the set of positive and negative energy states at 
this initial time. At this time the fluctuation operators may be expanded in this basis. The 
Fourier components of the fluctuation operators are thus written as 


Xjf(0) 




1 

vm 

i 


a fc(°) - 4(°). 
%(0)-4(0)] . 


J 


The number operators at any time t are 


N x,k(t) 

= 

N a ,k(t) 


Tral(0)admt) 

Tr p(0) 

r nr6t(°)6;(0)p(i) 

Tr p(0) 


b j ;(t) = U{t)b i: (0)U- 1 {t), 


Tra].(t)a £ (t)p(0) 
Tr p(0) 

Tr p(0) 


with U(t) the unitary time evolution operator. Following the arguments presented in refer¬ 
ence |TB[ we hnd that the creation and annihilation operators at time t are related to those 
at the initial time t — 0 by a Bogoliubov transformation. In terms of the dimensionless 
variables introduced above we hnd 


N\,q{ r ) — (2JF — 1) N Xj q(0) + (JF — 1) , 


F= l \U q {r)\ 2 


1 + 


\U q (r )\ 2 

wm q {r)\\ 


1 

+ 2 ’ 


N a ,q(r) = (2g - 1) N a>q (0) + {g - 1) , 


G = 7[ \V q (T) I 2 


1 + 


\V q (r )\ 2 
w 2 \V q (r)\ 2 


1 

+ 2 ’ 


where A^ x cr (0) are the occupation numbers at the initial time t = 0 and in the case given by 
the Bose-Einstein distribution functions, derivatives are with respect to the dimensionless 
variable r. Since in this section we will be working at zero temperature, only the last term 
will contribute to particle production. This term is recognized as the “induced” contribution. 
It can be seen from the initial conditions on the wave functions that the induced contribution 
vanishes at the initial time. 

We will compute numerically the number of particles produced in a correlation volume 


J\f b (t) = j d 3 kN k (t)/\m ^\ 3 . 


(3.7) 
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B. Inflaton Coupled to Fermions Only 


The fermionic Green’s functions are constructed from the solutions to the homogeneous 
Dirac equation in the presence of the background field. The main ingredients and treatment 


is similar to that of fermions in presence of an electric field studied by Kluger et. al. |23 
Writing the independent solutions as 

U^ 2 \x,t) = e i ^u£' 2 \t) , 

V ( 1 ’ 2 ) (f,t) = e~ its V k {h2 \t) , 

the mode functions obey 


*7o 


?7o 


dt 

d 

dt 


7 • k - m^(t) 
7 • k-m^(t) 


U‘i' 2 \t)= 0, 

vf- 2 , (i) = 0 . 


It is convenient to write the spinors as 


U, 


D’m = 


V7 2) (i) = 


- 7 • k + m^{t) 

d _ r 
+ 7 ' k + m ^\ t ) 


fk(t)u^’ 2) 

9k{t)v {1 ' 2) 


with u( 1,2 \ v^ 1,2 ^ the spinor eigenstates of 70 with eigenvalues +1, —1 respectively. The 
functions fk(t), gu{t) obey the second order equations 


dt 2 


+ k 2 + m 2 ,(f) 




fk(t) = 0 , 


dt' 2 


+ k 2 + m 2 ,(t ) + irn^it) 


9k(t ) = 0 . 


We now need to append initial conditions. We will consider the situation in which the 
system was in equilibrium at time t < 0 with the expectation value of scalar field being 0 ( 0 ) 
and 0(0) = 0. Thus the fermion mass is constant and given by m^( 0). We can now impose 
the condition that the modes fk(t), gkif) describe positive and negative frequency solutions 
for t < 0 and, normalizing the spinor solutions to the Dirac equation to unity, we impose 
the following initial conditions 


fk{t < 0) 
9k(t < 0) 


g—ikot 

\J2ko{ko + 777 ^( 0 )) 

gikot 


^j2k o (k o + m^(0)) ’ 


k 0 = 
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Equations (|3.8|) with these boundary conditions imply that 


9k(t) = . (3.8) 

The necessary ingredients for the zero temperature fermionic Green’s functions are the fol¬ 
lowing 


s;(t,t') = -i y. uirnK?), 

a=l,2 

siit,t') = i Y v°- k (t)v%t') . 

a=l,2 


In the standard Dirac representation for the 7 matrices, we hnd 

'1 + 7o 


s k (t, 0 = kit) flit') [W 4 (t) 7 o -7 -k + m^(t) 
Sk(t, t') = —ifl(t)f k (t') [Wfc(t)70 - 7 • * - "VW 


2 

1 - 7o 


Wfc(07o-7-^ + ^(0 
W fc (t')7o -7 -k-m^it') 


mft) = i 


Jk(t) 


fk(t) 


With these ingredients the non-equilibrium fermionic Green’s functions can be constructed 
as in section II. It is an important check that the equality Tr Sf + it, t) = Tr S k (£, t) = 
Tr Sf it, t) = Tr Sf(t, t) is fulhllcd where the trace is over Dirac indices. It is also an 
important property that 

(-ifk(t) + m i>ff)fl{t )) + m^(t)f k (tfj + k 2 flit)f k it) = 1 . (3.9) 

This is a consequence of the conservation of probability and may be checked explicitly. 

The evolution equation becomes 


'(pit) + ml (pit) + ^0 3 (t) + ^(pit) < X 2 ix,t) > -yTrS < (x,t]x,t) = 0 , 
o 2 


Tr S K (x, t; x,t) = 2 


d 3 k 

(2vr) 3 


l-2k 2 \f k {t)l< 


(3.10) 


The integral in ( |3.10|) is divergent. The divergence structure may be understood in two 
different ways: carrying out a WKB expansion of the modes, just as was done in the bosonic 
case (this is a tedious and lengthy exercise), or alternatively calculating the closed loop in the 
presence of the background held. We carried out both methods using an upper momentum 
cutoff in the integrals and found 


2 y 


d 3 k 


l-2k 2 \f k (t)f 


2n 2 


y A 2 m^(t) + -^[rdipit) + 2 m 3 (t)] In ^ + finite , 


(3.11) 


where k is an arbitrary renormalization scale which will be chosen again to be k = for 
numerical convenience. We will concentrate on the “chiral limit” in which the bare mass 
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term for the fermions vanishes and mp(f) = y(f>(t). In this limit the discrete symmetry 
0 —> —0 is explicit. The divergent terms in ( [3.1 ID are then identified as: mass renormaliza¬ 
tion (of 0) (first term), wave function renormalization (second term) and coupling constant 
renormalization (third term). From now on, all quantities will be renormalized, and we drop 
the subscript R to avoid cluttering. For simplicity in the numerical analysis and in order to 
isolate the fluctuations from the fermionic contribution from those of the scalar sector, we 
now neglect the contribution from the scalar loop. The scalar fluctuations had already been 
analyzed in the previous section and in [|l3j]. The final renormalized equation of motion is 
now 


A r d 3 k 

0(f) + m| 0(f) + -0 3 (f) - 2yJ — 
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An 2 - 


0(f) + 2 y 2 0 3 (f) 
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With the purpose of a numerical analysis of this equations, it proves convenient to introduce 
the following dimensionless of variables 


77 = r/0/m$, r = m$f , q=k/m$ >) g' = X R /6y 2 , g = y 2 /ir 2 , 

u q (r) = f k (t)^/uj 0 k (uj 0 k + y 2 (j) 2 ( 0)) , 


in terms of which the equations of motion are 


j^ + v + g'v 3 ~ 9Mr) = 0 , 


dr 2 
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£(r) = / q 2 dq 
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\Jq 2 + V 2 {0)(\Jq 2 + V 2 (0) + V(0))_ 
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r) + j (fj + 2 ry 3 ) In 
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rriR,<s> 
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sjq 2 + p 2 (0) 

V2 
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(3.12) 


Although the equations are finite when the ultraviolet cutoff is taken to infinity and can be 
written without the introduction of the cutoff by subtracting the integral with a lower limit 
cutoff given by the renormalization scale, numerically these integrals will have to be done 
by introducing an upper momentum cutoff anyways. Therefore we keep this UV cutoff in 
the equations. 


1. Particle Production 

The study of particle production is similar to the scalar case, with the only complication 
being the spinorial structure of the fermionic fields. At the initial time f = 0 the quantized 
fermion operator is 
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(3.13) 


2(20)= / y—Jj £ 

J V Z7V ) 0 = 1,2 

in terms of the mode functions determined above. The number of fermions and antifermions 
are defined as 


<Jv'(j)> = J 

r d 3 k ^ Tr 

1 (2^) 3 V 

[p(t)6i Q)t (0)^(0) 

r d 3 k ^ Tr 
' (2tt) 3 

p(°)4 Q) V)4 Q 4) 

Tr p(0) J 

Tr p(0) 

(N’(t)) = j 

r d 3 k y. Tr 

1 W 3 \ 

[ P (t)4 a)t (o)4 a) (o) 

r d 3 k Tr 

/ (2vr) 3 ^ 

P (o)4 a)t (t)4 a) (t) 

Tr p(0) 

Trp(O) 


The time dependent coefficients are obtained by projecting the time dependent spinor so¬ 
lutions onto the positive and negative energy solutions at t — 0, and are related to the 
coefficients at t — 0 via a Bogoliubov transformation 

t\t) = £ [rl;l + (t )b p . 
a 

= £ [«S fc+ («)r + . 

0 

The identities 

£ i <) i 2 +1 ?$k,-(t) l 2 = i, 

I 3)k,+ ( t ) I 2 + I ! 2= 1 

0 

ensure that the transformations preserve the anticommutation relations as they must. It 
is also a matter of algebra using the equations of motion for the mode functions to prove 
that the number of fermions minus antifermions is conserved for each k mode. In what 
follows, we only consider the number of fermions produced since fermions and antifermions 
are created in pairs. 

With some algebra, the Bogoliubov coefficients can be expressed as the function of f k (t) 
and flit). In terms of the dimensionless variables defined above we obtain the number of 
fermions within a correlation volume N*{t) = N* (r)/ 

^2 J Q 2 dqNj(T) 



x ^ + dir) - \jq 2 + ?/ 2 (0) - r/(0) u q (r) . (3.14) 
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IV. NUMERICAL ANALYSIS 


The numerical analysis was carried out with a fourth order Runge-Kutta method for the 
differential equations and a 5-point Bode rule integrator for the /c-integrals. The typical step 
size in time was 1 — 2.ICC 3 , and the typical step size in k was 10 -3 . The cutoff A/|m$| was 
varied between 75 and 200; we found no sensitivity to the cutoff in the range of parameters 
that we tested (see below). The code was very stable within the range of parameters tested. 


A. Inflaton Coupled to Scalars: Unbroken Symmetry Case 


Since the contribution of the one-loop quantum fluctuations of the inflaton <f> have already 
been studied previously []TT5[ and we want to study the contribution of the lighter field a, we 
will neglect the contribution from the self-interaction in the evolution equation. Thus we 
study the following equations obtained from equations (|3.3| , |3.G|) in terms of dimensionless 
variables 


V(r) +V(r) + r? 3 (r) + rj(r) J 1 ^ q 2 dq 


2 j lK ( r )| 2 ~ 1 
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+ %^ T ) [v 2 (r)-v 2 (0) 
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, % / N 

+ y w 


V q {r) = 0 , 


14(0) = 1 , 14 ( 0 ) = -iWg = -i 
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q 2 + 


m* 6 g 




2 + y^ 2 (°) • 


Notice that the factor (3/2) multiplying the logarithm in eq. (|3.6|) and missing from (|4.1|) 
arises from the renormalization of the <f>-scalar loop that is not taken into account in ( |4.1| ). 

Figures (l.a-c) show the unbroken symmetry case with the quantum fluctuations from 
the scalar loop of the a particles for the values y — 0; A/87r 2 = 0.2; g = A; rn a = 

0.2m$; 77(0) = 1.0; f/(0) = 0. Figure (l.a) shows r/(r) vs r, figure (l.b) shows A/’o-(t) vs r 
and figure (l.c) shows vs q for r = 120, similar graphs were obtained with snapshots 
at different (earlier) times. 

Figure (l.a) shows a very rapid, non-exponential damping within few oscillations of the 
expectation value and a saturation effect when the amplitude of the oscillation is rather 
small (about 0.1 in this case), the amplitude remains almost constant at the latest times 
tested. Figure (l.a) and figure (l.b) clearly show that the time scale for dissipation (from 
figure (l.a) is that for which the particle production mechanism is more efficient (figure 
(l.b)). Notice that the total number of particles produced rises on the same time scale as 
that of damping in figure (l.a) and eventually when the expectation value oscillates with 
(almost) constant amplitude the average number of particles produced remains constant. 
These figures clearly show that damping is a consequence of particle production. At times 
larger than about 40 mj 1 (for the initial values and couplings chosen) there is no appreciable 
damping. The amplitude is rather small and particle production has practically shut off. If 
we had used the classical evolution of the expectation value in the mode equations, particle 
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production would not shut off (parametric resonant amplification), and thus we clearly see 
the dramatic effects of the inclusion of the back reaction. 

In this unbroken symmetry case linear relaxation predicts a slow t~ 3 power law decay to 
an asymptotic finite amplitude because one particle decay is kinematically forbidden: the self 
energy contribution is a two-loop effect with one $ and two a particle cut and kinematically 
there is a one-particle pole below the three particle threshold. A slow power law linear 
relaxation asymptotically cannot be ruled out numerically because we have not continued 
the integration for longer times but clearly asymptotically the numerical result is compatible 
with linear relaxation. Figure (l.c) shows the distribution of particles created at the latest 
time r = 120 as a function of wave vector. Similar graphs were obtained with snapshots 
at different earlier times. The distribution is clearly non-thermal and skewed towards small 
momentum (in units of m$). These figures are qualitatively similar to those obtained in 
the self-interacting case in [|l3j. The asymptotic behavior is that of undamped oscillations 
of small amplitude. This is compatible with the result from the linear relaxation analysis 
because there is a one-particle pole below the three particle threshold resulting in undamped 
oscillations at large times. Linear relaxation predicts qualitatively the same results for the 
self-interacting scalar case and this case in the unbroken phase. The numerical results are 
consistent with this prediction. 

However, for large amplitudes non-linear relaxation via particle production is very effec¬ 
tive and dramatically different from linear relaxation. 


B. Inflaton Coupled to Scalars: Broken Symmetry Case 


In this case, we take m| = — |m||. 

As in the unbroken symmetry case, we only study the effect of the lighter a fluctuations. 
In the broken symmetry case, linear relaxation predicts an open decay channel for the 
inflaton, resulting in exponential relaxation for a long time and eventually relaxation with 
a power law. Therefore we should not expect a constant amplitude asymptotically but an 
amplitude that eventually should relax to zero. In this case the renormalized equations for 
evolution are 


v(t) ~ v{t) + h 3 0 ) + v{t) f'™*' q 2 dq 
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Figure (2.a-c) show rj(r) vs r, A 4(t) vs r and N q ^(j = 200) vs q respectively, for X/8n 2 = 
0.2; g/X = 0.05; rn a = 0.2 |m$|; g(0) = 0.6; fj( 0) = 0. Notice that the mass for the 
linearized perturbations of the held at the broken symmetry ground state is \/2 |m$| > 
2 m a . Therefore, for the values used in the numerical analysis, the two-particle decay channel 
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is open for linear relaxation. For these values of the parameters, linear relaxation predicts 
exponential decay with a time scale r re i ~ 300 (in the units used). Figure (2.a) shows very 
rapid non-exponential damping on time scales about six times shorter than that predicted 
by linear relaxation. The expectation value reaches very rapidly a small amplitude regime, 
once this happens its amplitude relaxes very slowly. Within our computing time limitations 
we could not confirm that there is exponential relaxation in the small amplitude regime (for 
r > 100) but clearly there is a striking difference with the unbroken symmetry case. The 
influence of open channels is evident, however in the non-linear regime relaxation is clearly 
not exponential but extremely fast. Although we cannot confirm the exponential (or power 
law) relaxation numerically in the small amplitude regime, the amplitude at long times seems 
to relax to the expected value, shifted slightly from the minimum of the tree level potential 
at rj — 1. This is as expected from the fact that there are quantum fluctuations. Figure 
(2.b) shows that particle production occurs during the time scale for which dissipation is 
most effective, giving direct proof that dissipation is a consequence of particle production. 
Asymptotically, when the amplitude of the expectation value is small, particle production 
shuts off. We point out again that this is a consequence of the back-reaction in the evolution 
equations. Without this back-reaction, as argued above, particle production would continue 
without indefinitely. Figure (2.c) shows that the distribution of produced particles is very 
far from thermal and concentrated at low momentum modes k < |m$|. This distribution is 
qualitatively similar to that in the unbroken symmetry case, and points out that the excited 
state obtained asymptotically is far from thermal. 


C. Inflaton Coupled to Fermions Only: Unbroken Symmetry Case 


Here we treat the case y ^ 0; g = 0. 

The renormalized evolution and particle production equations in this case are given by 
eq. ( |3.12| ) and eq. (|3.14|) , respectively, in terms of dimensionless quantities. Figures (3.a-c) 
show 7](t) vs. r, vs r and J\fj (r = 200) respectively for the values of the parameters 

m^ = 0; y 2 /tt 2 = 0.5; A/6 y 2 = 1.0; 77 ( 0 ) = 0.6; ?)(0) = 0. One observes from these 
figures that after a rather brief period of initial damping of just a few oscillations of the scalar 
field, the dissipative mechanism shuts-off. Figure (3.b) shows that during this time scale 
fermion-antifermion pairs are being produced but then the number of produced particles 
saturates and oscillates with a small amplitude. Figure (3.c) shows that the distribution of 
fermions produced is peaked at very low momentum (k < m$) and with a maximum value 
of 2, which is the total number of degrees of freedom per k -wave vector. These numerical 
results expose the physics of Pauli blocking very clearly; the available low momentum modes 
are occupied and no more fermion-antifermion pairs can be produced. Pauli blocking shuts 
off particle production and dissipation very early on. We have obtained snapshots of the 
particle number as a function of momentum for different times, and they all present the 
same picture. 

This result is markedly different from the prediction of linear relaxation. At zero temper¬ 
ature, eq. (|2.48|) for linear relaxation predicts exponential damping with a (dimensionless) 
time scale T re i ~ 10.6 for the values of the parameters chosen above. The difference between 
the non-linear evolution and that predicted by linear relaxation is explained by the Pauli 
blocking phenomenon. Very early in the evolution, the low momentum available fermionic 
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states were filled with produced fermions. Once these states have filled, damping and par¬ 
ticle production shuts off. This Pauli blocking effect is explicit in eq. (|2.48|) but there it 
appears from finite temperature effects. In the non-linear relaxation case, we began at zero 
temperature, but an excited state quickly ensues because of particle production. 

This analysis reveals that for large amplitudes of the scalar field, fermions will be rather 
ineffective in dissipation and damping because of Pauli blocking even at zero temperature. 
The time scales for dissipation obtained from the fermion self-energies, which apply to the 
case of linear relaxation are completely unrelated to the time scales for non-linear dissipative 
processes, even asymptotically, because the fermionic states are Pauli blocked. 


V. CONCLUSIONS AND IMPLICATIONS 

The reheating and thermalization processes in inflationary universe models occur very 
far from equilibrium and must be studied in their full complexity, eventually numerically. 
The methods developed within real-time non-equilibrium field theory allow a consistent 
treatment that we used to obtain the equations of evolution for the expectation value of the 
inflaton field. These equations are non-perturbative and take into account the back reaction 
effect of both the self couplings of the inflaton, as well as the couplings to lighter scalars and 
fermions. 

We have examined these evolution equations, both in the linear regime, in which the 
amplitude of the field is small and in the opposite extreme, the non-linear regime, which is 
intrinsically non-perturbative. 

In the linear regime, we have seen that if the inflaton mass is above the two particle 
threshold, there is some damping behavior for a period of time, and indeed, this damping 
is governed by the total width of the inflaton. However, at late times the behavior is 
dominated by power law decay rather than exponential damping. Furthermore, we have 
noted the important distinction between the time scale for the relaxation of the expectation 
value of the inflaton as opposed to thermalization of the produced particles. It is clear that 
as treated in this work, relaxation can occur well before thermalization, i.e. interactions that 
drive the particle distributions towards a thermal one, has a chance to become relevant. 

The more interesting case, and the one which would most likely be relevant to the 
reheating problem in models such as chaotic inflation j5j, is that of when the inflaton is in 
the non-linear regime. Unlike the elementary approach to reheating, in which the “decay” of 
the coherent inflaton oscillations occured a time tdecay ~ T _1 , where T is the inflaton decay 
width, after the end of the slow-roll regime, particle production in the non-linear regime 
begins at very early times and then shuts off. Indeed, by the time the inflaton expectation 
value reaches its asymptotic state of oscillations around its minimum, the period of particle 
production is essentially over. This will be the stage of thermalization via collisions, however 
we find that the spectrum of produced particles is extremely non-thermal, and collisional 
relaxation towards an equilibrium thermal state may take a long time. 

An inflationary model that cannot reheat the universe to at least nucleosynthesis temper¬ 
atures is wrong. Thus it is clearly important to be able to compute the reheating temperature 
in such a model. What our analysis here shows is that this computation is of necessity much 
more complicated than has previously been thought but perhaps more importantly that 
there are different time scales. 
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The first stage which occurs rather fast is that of particle production and damping of 
oscillations as a result of induced amplification. The second stage, that begins typically 
when the amplitude of the inflaton is rather small and it oscillates with almost constant 
amplitude at the bottom of the potential is that of one-particle decay and of thermalization 
via collisional relaxation. 

The thermalization time will have to be studied setting up a (quantum) Boltzmann 
equation using the distributions of produced particles at the end of the induced amplification 
stage as input for the Boltzmann evolution. Since these initial distributions are very far 
from equilibrium, it may take many collisions to relax to a thermal equilibrium state. If 
the couplings are very small (and this will clearly depend on the models) the thermalization 
time may be very large. 

Following the particle numbers to the point at which they become thermal would then 
allow us to pick off the final reheating temperature. In the weak coupling regime, as men¬ 
tioned above, the thermalization rate may be smaller than the expansion rate, allowing 
for significant redshifting of the total energy density, leading to a low reheat temperature, 
perhaps of the order of the inflaton mass or less. 

The particular case of large thermalization rates compared to the expansion rate, allows 
for a “quick” estimate of the reheating temperature. In this case we can assume that 
thermalization occurs without red-shifting of the energy which is then conserved during the 
thermalization time. By taking the energy density at the beginning of this stage to be aT A 
with a the Stephan Boltzmann accounting for the particle statistics one can then obtain an 
estimate of the reheating temperature, but the justification for this should ultimately arise 
from a deeper understanding of the time scales involved. 

Our main point in all of this, though, is that a more detailed calculation will have to be 
done in order to extract the reheating temperature in a more reliable way. However, what 
we have done is to clarify what the relevant particle production mechanisms are in these 
theories and to provide a consistent and implcmentable framework of calculation. We are 
now extending these studies to FRW cosmologies. 

We can compare our work to other recent work in this subject. Brandenberger et al [12 
have performed a calculation of particle production for an oscillating inflaton held. However, 
they have used the inflaton held as a background which does not respond to the produced 
particles. In other words they neglect the back reaction, which as we have seen, allows for 
the production process to shut oh, and determines the time scales for particle production. 
Kofrnan, Linde and Starobinsky El have also performed such a calculation, and where we 
have common values of the parameters, we agree. However, they do not follow the produced 
particles through the thermalization period, so that we are somewhat uncertain about how 
they draw their conclusions about the final reheating temperature. We emphasized in this 
study that thermalization and particle production are fundamentally different processes and 
in the non-linear regime likely to occur on widely different time scales. 

One area of current interest in which this work may have some implications is that 
of the so-called “post-modern Polonyi problem” concerning flat directions for some of the 
moduli fields in string theories |17| . These are holds with pertubatively hat directions whose 
degeneracy is lifted by non-perturbative effects. Their masses then become of order the 
weak scale and their couplings to normal matter are gravitational. These properties allow 
the energy density in these helds to dominate that of the radiation in the universe until 
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times well after nucleosynthesis. 

This may not necessarily be true given our analysis. Since particle production occurs 
not just during oscillations around the minimum now but throughout the evolution of the 
the moduli held, it is not impossible that some of this energy density could be transferred to 
lighter particles, prior to nucleosynthesis, perhaps averting this potential catastrophe. We 
are currently looking into this possibility. 

To reiterate then: induced amplification during the evolution of the inhaton held allows 
for a very different mechanism for particle production at the end of reheating than has 
been used in the past for inhationary models. In the non-linear regime, this is by far the 
most important such process, and will have significant consequences for inhationary universe 
models. 
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APPENDIX A: A PEDAGOGICAL EXERCISE 


In this appendix we show explicitly how to implement the tadpole method for the case of 
the inhaton coupled to a lighter scalar held via a simple trilinear coupling. We will not worry 
about renormalization issues here, though they were, of course, dealt with in the text. The 
formalism for non-equilibrium quantum held theory has already been described several times 
in the literature [^6|. A path integral representation involves an integration along a path in 
the complex time plane with forward, backward branches and if the initial density matrix 
was that of an equilibrium system at an (initial) temperature T also an imaginary time 
branch. For real time correlation functions the imaginary time branch does not contribute 
and only determines the boundary conditions on the Green’s functions. 

We use the tadpole method @,|I3|^] to study the time evolution of the expectation value 
of the inhaton. 


Mx,i) = Tl ' _ Tr [4>-(a?)p(t)] 


(Al) 


Trp(O) Trp(O) 

where are the helds defined on the forward and backward branches respectively. We set 


$ ± (x, t) = <t>(x, t ) + x ± (x, t ) , (A2) 

where y 1 * 1 are held fluctuation operators defined along the respective branches. 

The effective evolution equation for the background held 4>(x, t ) follows from the condi¬ 
tion. 
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< X ± (x,t) >= 0 . 


(A3) 


Treating the linear and non-linear terms in as interactions and imposing the condition 
(|A3j) consistently in a perturbative or loop expansion, one obtains expressions of the form 
(here we quote the equation obtained from < x + (T, t) >= 0 ) 


dx'dt' [< y+(f, t)x + {x ', t') > 0 ++ {x', t')+ < x + (x, t)x~(x', t’) > 0 + ~(x', t') 1 = 0 (A4) 


and similarly for < x~{x,t) >= 0. The 0 ±i: are in general integro-differential operators 
acting on the background held. 

Because the Green’s functions < x + (x, t)x + (x', t') >, etc. are all independent one obtains 
the equations of motion in the form 

0 ++ (x',t') = 0 , 0 + ~(x' ) t')= 0, 0~ + (x',t') = 0, 0~~(x',t') = 0 . (A5) 


It is a consequence of the properties of the non-equilibrium Green’s function (see below) 
and ultimately a consequence of unitarity that all the integro-differential operators O are 
the same. 

Consider the Lagrangian density 

£(<h,a) = £ 0 (<L) +£ 0 (o-) + g(t) $ cr 2 , (A6) 


with the £o being the free held Lagrangian density (with respective mass terms) and 
have allowed the coupling to depend on time. The non-equilibrium path integral requires 
£( < h + ,a + ) — £(<I) - ,<t - ). In the tadpole method we write <h ± (x, t) = x ± (x,t) + cp(t) and 
identify <p{t) as the (non-equilibrium) expectation value of the held $ . This identification 
then requires that < x(x,t) >= 0 where the expectation value is in the non-equilibrium 
density matrix with the path integral representation along the contour in complex time. 
After this shift, the action reads: 


L 


d 3 xdt {£ 0 (X + ) + £ 0 (<t + ) + X + 

(x + ,o’ + -*■ • 


-0 - + g(t)4>(t)(a + ) 2 + g{t)x 


+ (a + ? 


The linear term in is included as a perturbation. The hrst contribution to the equation 
of motion is obtained from this linear term, from the condition < x + (x, t ) >= 0 one obtains 
to this order 


J dt' {< x + (t)x + (t') > 00 [-0 - m|0 


- < x + (t)x (0 > 0) \-4> - m 


= o 


(AT) 


where we have suppressed the spatial arguments. Because the correlation functions 

< X + (t)x + (0) >, < X + (t)x~(t') > are independent, one obtains the tree level equa¬ 

tions of motion. It is straightforward to see that the same is obtained by imposing 

< x~{x,t) >= 0. In an amplitude expansion (an expansion in powers of 0(t)) the one 
loop correction to the equation of motion is obtained by expanding (the exponential of) 

g0(t)(a + ) 2 + gx + ( a+ ) 2 — (H -> —) to hrst and second order. The hrst order gives a tadpole 

contribution, the second order needs one vertex with y, the other with 0. One obtains 
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dt' < X + (f)y + (f') > {(*) (-0 - m|0) + (ifl'(O) < (o’ + (0) 2 > + 

' dt"(ig(t")) 2 [< ( a + ( i ')) V + (*")) 2 > - < (< T + («')) VW ")) 2 >] •#>(«")} - 
df' < X + (0)x~CO > {(0 (-0 - m|0) + (iff(t')) < (O 2 (0 > - 
’dt"(-ig(0) 2 \< (a~(t')) 2 (a-(t")) 2 > 0(f")- < {a~ (f'))V + (f)) 2 >1 0(0 


= 0 


The expectation values are computed using Wick’s theorem and using the free-held Green’s 
functions of section II. The tadpole (time independent) is absorbed in a shift of the expec¬ 
tation value. The coefficient of < x + X + >> < X + X~ > must vanish independently because 

these Green’s functions are independent and must vanish at all times. From the Green’s 
functions (and more generally from the formal time contour integral) it is seen that the 
equations obtained are identical. If the expectation value is translational invariance, the 
spatial integrals set the momentum transfer to zero. For example using the zero tempera¬ 
ture Green’s functions of section I, one finds that the term that has the non-local (in time) 
correlation functions (last term) becomes 


dt'\ig{t")) 2 < 0r+(f'))V+(f")) 2 > - < ^ + {t')) 2 {a-{t")f > 0(0 = 


2 / df"(i<?(f")) 2 0(f")©(f' - t" 


d 3 k f _ 2 ^sin [2w k (t' - f")] 


(2vr 


*4 


For the resummed one-loop approximation the term ^0(t)(cr ± ) 2 is absorbed in the mass 
term for the a held and only the ycr 2 terms are considered as perturbation. In this case 
the 1-loop contribution is obtained from the tadpole term < (cr ± (t)) 2 >, which is now time 
dependent and obtained from the non-equilibrium Green’s functions constructed with the 
mode functions for the time-dependent mass as in section IV. One now finds the evolution 
equations 


0(f) + ml 0(f) + g{t) < v 2 {x, t) >= 0 , 

Ai.t) >=-> j =I " m) ' 2 


< 


(2vr)3 2w k 


coth 


(3w k 


d 2 


^ + k 2 + m 2 (t ) | 14(f) = 0 , 


14(0) = 1 , 14(0) = -iw k = -i\Jk 2 + m 2 (0) , 

m 2 (f) = m a + g(t)4>(t) . 
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Figure Captions: 

Fig.(l.a) Scalars: unbroken symmetry t](t) vs t for the values of the parameters 
y = 0; A/87T 2 = 0.2; g = A = 1; m a = 0.2 m^; 77 ( 0 ) = 1.0; r/(0) = 0. 

Fig.(l.b): A f a {r) vs. t for the same value of the parameters as figure (l.a). 

Fig.(l.c) : -V 9i(T (r = 120) vs. q for the same values as in figure (l.a). 

Fig.(2.a) Scalars: broken symmetry ?/(r) vs r for the values of the parameters y = 
0; A/8tt 2 = 0.2; g = A = 0.05; m a = 0.2 |m^|; 77 ( 0 ) = 0.6; t)(0) = 0. 

Fig.(2.b): A f a (r) vs. r for the same value of the parameters as figure (2.a). 

Fig.(2.c) : A fq, a {j = 200) vs. q for the same values as in figure (2.a). 

Fig. (3. a) Fermions: unbroken symmetry t](t) vs r for the values of the parameters 
9 = 0; y 2 /ir 2 = 0.5; A/6 y 2 = 1; = 0; 7/(0) = 1.0; t)(0) = 0. 

Fig.(3.b): AC(r) vs. r for the same value of the parameters as figure (3.a). 

Fig.(3.c): A/V(r = 200) vs. q for the same values as in figure (3.a). 
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